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Abstract 

I observe a rough volume-dependent scaling of the low eigenvalues of a chiral Dirac operator in 
lattice studies of SU(3) lattice gauge theory with two flavors of color sextet fermions, in its weak- 
coupling phase. The mean value of the ith eigenvalue scales with the simulation volume V = L'^ as 
{Xi)L'^ ~ Q, where Q is a volume-independent constant. The exponent p is about 1.4. A possible 
explanation for this phenomenon is that p is the leading relevant exponent associated with the 
fermion mass dependence of correlation functions in a theory whose zero-mass limit is conformal. 



* Electronic address: |degr and@pizero.colorado.edu 




FIG. 1: Average values of 8 lowest eigenvalues of the sextet-representation valence overlap operator, 
vs where the lattice volume V is defined to be y = L^. The actual volumes, moving from 
left to right across the graph, are 16^, 16^ x 8, 12^, 12'^ x 8 and 12'^ x 6. The four lines are scaling 
curves {\i)LP = constant, for p = 4, 3, 2, 1, from the top of the figure down. 



I. INTRODUCTION 



Recently, various research groups [H, S, i, i, H, 0, 0, S, i, E Ell, E, [3, 0, H [13 have 
begun to use lattice methods to study field theories which might be candidates for strongly 
coupled beyond - Standard Model phenomena [l7|. These models typically involve gauge 
fields and a large number of fermion degrees of freedom, either many flavors of fundamental 
representation fermions (where the discussion goes back to Refs. [l8,[l3) or a smaller number 
of flavors of higher- dimensional representation fermions (strongly emphasized by Refs. (2^, 
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I have been part of a collaboration studying SU(3) lattice gauge theory with two flavors of 
color sextet fermions [sl, 11 1. This theory has a weak-coupling phase which is deconfined and 
chirally restored. Recently, I observed a peculiar volume-dependent scaling of the eigenvalues 
of a valence quark Dirac operator, in the weak coupling phase. Fig. [T] shows this result. 
This is very different from the spectrum in simulations with fundamental representation 
fermions. A presentation of the dependence, a comparison with simulations of two flavors 
of fundamental representation fermions (which do not show this scaling dependence), and a 
possible explanation, are the subjects of this paper. 

I begin with some (possibly) relevant background: The usual description of renormaliza- 
tion for a gauge theory coupled to massless fermions generally concentrates on the behavior 
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of the running gauge coupling, and makes a distinction about how the gauge couphng flows 
under rescaling to the infrared: it could flow to strong coupling in the case of a confining 
theory, or to zero, for a trivial theory, or to a fixed point at some nonzero value, g^* . The 
latter case is referred to as an infrared-attractive fixed point (IRFP) theory. This is a phase 
with no confinement, no chiral symmetry breaking, and algebraic decay of correlation func- 
tions (no particles). Evidence has been presented that several such theories exist: SU{'i) 
gauge theory with Nf = 12 flavors of fundamental fermions 3, [13, SU{3) gauge theory with 
Nf = 2 flavors of sextet fermions (the subject of this paper) |5| and SU{2) gauge theory 
with Nf = 2 flavors of adjoint fermions (l5|. 

While for confining theories the gauge coupling is relevant (the Gaussian fixed point g = 
is unstable), that is not the case for IRFP theories. The gauge coupling (more precisely the 
distance of the bare gauge coupling from its fixed point value {g'^ — g"^*)) is an irrelevant 
coupling. The critical surface encompasses a wide range of values of bare gauge couplings, 
plus the values of any other irrelevant operators (in a lattice theory, these are the usual lattice 
artifacts). The relevant coupling of an IRFP theory is the fermion mass, which must be fine- 
tuned at the UV scale to reach the critical surface (i.e. = 0). How a nonzero fermion 
mass drives the lattice theory off the critical surface determines the asymptotic behavior of 
(nearly) all correlation functions. The situation is completely equivalent to that of an order- 
disorder transition in a magnetic system. The only difference is that the relevant direction is 
parameterized by the quark mass rrig, rather than the reduced temperature t = {T — Tc)/Tc 
of the magnet. (A closer analogy is to a system which has been fine tuned to its Curie point 
and then placed in an external magnetic field. The external field breaks the underlying 
global symmetry just as a quark mass explicitly breaks chiral symmetry.) [1^ 

The framework to describe the physics of these systems is also standard. Tuning the 
mass to zero causes the correlation length to diverge algebraically, 

i_ 

^ ~ mg . (1) 

As m is tuned to zero, the singular part of the free energy scales as 

/,(m,) = m^/y-{A, + A2m[y-^\/y-). (2) 

where D is the system's dimensionality (here D = 4) Ai and A2 are non-universal constants 
and Hi is the biggest non- leading exponent. This is most likely the exponent i/g of the gauge 
coupling, g'^ — g'^*, which can be determined from the beta function as measured in (for 
example) Schrodinger functional simulations at mg = 0. 

How can one measure The most direct way is through the correlation length ^, 

through Eq. [H In a finite simulation volume, ^ will obey Eq. [1] until it grows to be the order 
of the size of the lattice, and then will saturate. To observe Eq. [1] probably requires large 
volumes. 



Unitarity bounds for conformal field theories [2J] constrain the critical exponent of the 
condensate, = S, 



m,) = ^ ~ (3) 



to lie in the range 3 > 7 = D/ym — 1 > 1. (Note that for D = 4, 7 is defined to be 3 — 7^ 
where 7^ is the anomalous dimension for the fermion mass.) For D = 4 this means that 
the allowed range of ym is between 1 (free field behavior, all exponents except yg given by 
dimensional analysis) and 2. This make a measurement of y^n through the mass dependence 
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of the condensate problematic. The condensate has a UV- sensitive piece proportional to 
the quark mass 25|, and so {il)^)) ~ Ciniq + C2iTi^ + . . . . This masks the non-analytic 
behavior. This is quite similar to the situation in finite temperature QCD, precisely at 



Tc, where [26 



S(a, TUq, T) ~ Cimq/a? + c^m}^^ + analytic. (4) 

It is different from QCD in that for the IRFP theory there are no Goldstone bosons (which 
contribute their own non-analytic piece to the condensate, below Tc). It is also different 
from QCD in that while in QCD, Eq. Hlapplies only at Tc, in an IRFP theory Eq. Hlgives the 
behavior of the condensate throughout the basin of attraction of the IRFP. Worryingly, the 
expected result from D = 3 0(4) universality, 1/5 = 1/0.56 = 1.78, implies that the QCD 
number is smaller than the smallest unitarity bound, meaning that the QCD susceptibility 
would then be more singular than in a conformal theory. However, the QCD susceptibility 



does not seem to have been measured yet [27 



The way I choose to attack through the condensate involves the Banks-Casher 
relation [2^, which relates the condensate to the density of eigenvalues A of the Dirac oper- 
ator p(A). At nonzero mass it is 



Iraq 



- - - / "i^)''^^- (5) 



If the theory is conformal, and S(mq) ~ m?^ then p(A) also scales as A'''. The search for 



the effect I saw was motivated by a finite-size scaling argument [29j which relates the scaling 
for the density p to the scaling of the value of individual eigenvalues. If we consider the 
average value of the ith eigenvalue of the Dirac operator in a box of volume V = L^, and if 
p{X) ~ A"^, then we expect 

(A.) - (tY (6) 



where the exponent is 

P=T^- (7) 

1 + 7 

(A quick derivation: p ~ A''' means that eigenvalues are uniformly distributed in an = 7+I 
dimensional space of volume V = , 

N 

^ = \C£n^Y'' r^, = l,2,...i? (8) 

i=l 

SO an eigenvalue scales as Aj ~ 1/R = {1/V)'^ = {1/V)^. Now suppose we are in D 
physical dimensions; in a box of volume V, there are V = modes, from which Eq. [7] is 
obtained. One example of this formula is free field theory: 7 = D — 1 and p = 1. Another is 
the case of chiral symmetry breaking encoded in the usual formulas of its Random Matrix 
Theory analog: 7 = so p(A) ^ po a constant, and p = D. This is an eigenvalue spectrum 
which depends on the dimensionless product XHV or (Aj) ~ ^/V-) 

For the case of an IRFP theory, p is equal to ym, the leading exponent. A way to determine 
ym is through the volume dependence of eigenvalues of the Dirac operator. 

There are several problems associated with carrying out this proposal. The first is nu- 
merical. Lattice discretizations of the Dirac operator which are easy to implement (allowing 
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simulations on large volume) break the full SU{Nf) ® SU{Nf) chiral symmetry of the con- 
tinuum Dirac operator. This complicates the analysis. For example, for the fermions in 
my simulations (Wilson- type fermions), explicit chiral symmetry breaking in the action in- 
troduces an additive shift to the eigenvalues. For a cleaner test of Eq. El I will use partial 
quenching: I will take configurations generated with lattice fermions which do not have 
exact chiral symmetry, and measure the Dirac spectrum using an implementation of lattice 
fermions with exact chiral symmetry (overlap fermions). This raises another issue: is what 
the valence fermion sees a faithful realization of what is happening in the equilibrium dis- 
tribution of real dynamical variables? I think that for a first study, what I am going to do 
is adequate. 

The next problem involves interpreting the results. For the case of a system which 
exhibits chiral symmetry breaking, there is a tight theoretical description of the behavior of 
the lowest eigenvalues of the Dirac operator, which allows one to relate delicate features of 
the spectrum to the low energy constants of the theory (the condensate, the pseudoscalar 
decay constant, and possibly others). This description is based on Random Matrix Theory 
(RMT). It is quite straightforward to use Random Matrix Theory predictions to tell whether 
a system exhibits chiral symmetry breaking. Conversely, if a system is in a chirally-restored 
phase, there are no longer RMT predictions to compare results against. 

I already know that our target theory has a weak-coupling phase which is deconfined and 
chirally-restored. The restoration of chiral symmetry is observed through regularities in the 
spectrum of screening masses as well as the behavior of the pseudoscalar decay constant as 
a function of quark mass. I can only continue while using the simplest properties of the 
eigenvalues, namely their scaling with system size. 

Now I return to a fuller description of Fig. [H This is a plot of the average value of the ith 
eigenvalue of the valence overlap operator versus simulation volume, plotted using L 
The data are all collected at the same values of the bare parameters, in the deconfined phase 
and at a quark mass where the lattice volume causes the correlation length to saturate. The 
actual volumes, moving from left to right across the graph, are 16^, 16^ x 8, 12^, 12^ x 8 
and 12^ x 6. Clearly, the spectrum does not appear to depend on much beside the volume, 
and it seems to show scaling, (Aj) ~ L"^. 

The four lines are scaling curves (Aj)L^ = constant, for p = 4, 3, 2, 1, from the top of 
the figure down, p = 4 would be the scaling curve if chiral symmetry were broken, p = 1 is 
free- field behavior. The data appear to lie in between - in fact, the best fit value is about 
p = 1.4. 

Svetitsky, Shamir and I have already observed (but with a different lattice discretization) 
that ^^/(S) N f = 2 sextet QCD shows evidence for an IRFP [5|. If that observation survives 
future tests, the exponent p is indeed ym, the leading relevant exponent. If that observation 
turns out to be false, the volume scaling of eigenvalues is still something striking. It is 
different from what is seen in the deconfined phase SU (3) Nf = 2 fundamental QCD (as we 
will see, below). 

The data in Fig. [T] all come from simulations at one set of bare parameters. However, if 
the simulations were in fact done in the basin of attraction of the IRFP, irrelevancy of the 
gauge coupling means that different gauge couplings merely correspond to lattice actions 
with different amounts of scaling violations. 

Finally, the idea behind this calculation was first described by the authors of Ref. j^. 
They proposed doing simulations with dynamical chiral fermions (overlap fermions) so that 
the Dirac operator whose eigenvalues are measured is the same one that appeared in the 
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action. They performed these simulations, also in the deconfined phase of SU{3) Nf = 2 
sextet QCD. They identified that the eigenvalue spectrum was inconsistent with expectations 
for the behavior of a system with chiral symmetry breaking, but only studied one volume 
(6"^ lattices). 

The paper proceeds as follows. In Sec. [II] I describe details of the simulations. In Sec. HII] 
I provide some background: I examine whether lattice data show evidence that the gauge 
coupling is irrelevant. I make a direct attack on the correlation length exponent using Eq. [TJ I 
find numbers also in the range 1.5-1.6. I then make some comparisons of Nf = 2 fundamental 
and sextet QCDs in their deconfined phases. These systems are quite different. In Sec. IIVI 
I provide background for Fig. [1] and describe my attempts to pin down the exponent p. I 
conclude with some speculations and (rather obvious) suggestions for follow-up work. 



II. NUMERICAL TECHNIQUES 

I performed simulations on a system with SU{3) gauge fields and two fiavors of dynamical 
fermions in the symmetric (sextet) representation of the color gauge group. The lattice 
action is defined by the single-plaquette gauge action and a Wilson fermion action with 
added clover term sO]- The fermion action employs the differentiable hypercubic smeared 



link of Ref. [sit], from which the symmetric- representation gauge connection for the fermion 
operator is constructed. No tadpole-improvement is used and the clover coefficient is set to 
its tree-level value. The bare parameters which are inputs to the simulation are the gauge 
coupling f] = Q/g'^, the fermion hopping parameter k. The inte grat ion is done with one 



additional heavy pseudo- fermion field as suggested by Hasenbusch j32|, multiple time scales 
(33| . and a second-order Omelyan integrator |34l |. 

The routines for simulating sextet-representation fermions were developed with (and 
mostly by) B. Svetitsky and Y. Shamir. The dynamical fermion algorithm was adapted 
from a program written by A. Hasenfratz, R. Hoffmann and S. Schaefer|35]. All computer 
code is based on the publicly available code of the MILC collaboration [36] . 

Simulation volumes range from 12^ x 6 to 16^ sites, and typical data sets range from a 
few hundred to a thousand trajectories. I recorded lattices every five trajectories (of unit 
simulation time trajectory length) and collected 40-80 lattices per parameter set for the 
calculation of screening lengths and overlap eigenvalues. 

The trick of combining periodic and anti-periodic boundary conditions for valence quarks 
37 , 38, [3^, 13] is used in spectroscopy or screening mass measurements. 



Throughout this work, instead of quoting k, I will use the the Axial Ward Identity (AWI) 
quark mass, defined through 

dt J2 (M^^ t)Xm = 2m, J2 (^(^' ^)^(O)) ■ (9) 

X X 

where = ?/'7o75'?/', P = iplnip, and X is any source. (I temporarily drop factors of the 
lattice spacing in the derivations.) The derivative is taken to be the naive difference operator 
{di,f{x) = {f{x + fia) - f{x - fia))/{2a)). 

The valence Dirac operator whose eigenvalues are the subject of this study is the overlap 
operator (4ll . l4p| . Details of the particular implementation of the action are described 
in Refs. |43l . |4J, |45|, |46|, |47|. The only new ingredient is the application to symmetric- 
representation fermions, using the same combination of hypercubic link and projection as 
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for the dynamical fermions. Eigenvalues of the squared Hermitian Dirac operator D^D are 
computed using the "Primme" package of McCombs and Stathopoulos|48| and split apart 
in the usual way. 

There are potential problems with this analysis. The first one involves the index theorem, 
relating the winding number of the gauge field k to the number of Dirac fermion zero modes, 

index = 2T{R)k. (10) 

T{R) is the Dynkin index of the representation R, 1/2 for fundamental representation 
fermions, {N + 2)/2 for sextet fermions in the color group SU{N), and so on. Thus we 
expect to see multiples of 5 zero modes for adjoint overlap fermions in our SU{3) case. 

However, ten years ago Heller, Edwards, and Narayanan discovered [i^ that the index the- 
orem applied to adjoint overlap fermions in background SU{2) gauge configurations failed: 
here 2T{R) = 4 and they saw configurations with zero modes which were not multiples 
of four. Similar results were more recently reported by Garcia Perez, Gonzalez-Arroyo 
and SastreSOj. Nogradi reported similar behavior from simulations with SU{3) sextet 
fermions 5l|] and in simulations where the bare gauge coupling is large, I have seen con- 
figurations whose zero mode content was a not a multiple of 5. 

It is unknown whether this is a disease, just a particular failure of the overlap action to 
capture topology when the gauge configuration is rough, or real physics, something which 
persists in a continuum limit. I am going to ignore it. The gauge configurations at the 
parameter value used to construct Fig. [1] were smooth enough that all of the lattices I 
collected for Fig. [T] had Q = 0. 



The second potential problem concerns possible phase structure in the theory. In Ref. [11 
we observed that in the deconfined phase of sextet QCD, the Polyakov loop ordered in 
one of the negative real directions, roughly along one of the complex elements of Z{3) 
(Re(TrP(a:)) < 0, Im(TrP(x)) ^ 0). This is quite different behavior from fundamental 
QCD, where the Polyakov loop orders positively (Re(TrP(a:)) > 0, Im(TrP(a;)) = 0). It is 
likely that there is a complicated phase structure as a function of the bare parameters, with 
different vacua favored in different regions of the bare parameter space jH^] . Detailed studies 
of the phase structure (at smaller volumes) are currently being performed by Machtey, 
Shamir, and Svetitsky [53|. For the time being, I believe, based on my simulations, that 
the equilibrium small quark mass region is Re (TrP(x)) < 0. I checked that all the data in 
Fig. [T] remained in the Re (TrP(x)) < 0, Im (Tr P(x)) ^ 0, phase during their time histories. 



III. SOME SPECTROSCOPIC CHECKS 

Before analyzing eigenvalues, I pause to consider some potential complications. 



A. Does the gauge coupling appear to be irrelevant? 

If the gauge coupling is irrelevant, observables collected at different values of the bare 
coupling will show small differences, qualitatively similar to the different sizes of scaling 
violations seen in ordinary QCD simulations when different lattice actions are used. 

We jll| saw this behavior in our earlier work, but did not describe it as such. It is easy to 
overlay screening masses taken in the deconfined phase at different values of j3 (at identical 
lattice volumes) and observe that they coincide. I now show a variation of that comparison. 
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FIG. 2: Map of the bare coupling constant plane relevant to our sextet simulations. The solid 
line is the line of zero quark mass, k = Kc- The crosses show the location of the confinement- 
deconfinement crossover at Nt = 6. Octagons show the data collection points for the next figures. 

To begin, the phase diagram of the present system is shown in Fig. [2J The hne is the 
location where the AWI quark mass vanishes, k = k^. The crosses show the location of the 
Nt = & deconfinement transition. To the left of the crosses, the theory confines and (to the 
extent we can discuss this for massive quarks) chiral symmetry is broken; to the right, the 
theory is deconfined and chiral symmetry is restored. Octagons show the data collection 
points (all on 12^ x 6 volumes) for the next figures. 1 have checked that all data collected 
at larger volumes (all at /? = 5.2) is deconfined, but have not accurately mapped the phase 
boundary. 

1 define the correlation length ^ to be the inverse of the screening mass in the pseudoscalar 
channel and plot ^ vs l/ruq, the inverse AWI fermion mass, in Fig. [HJ (In the deconfined 
phase, the pseudoscalar, vector, axial vector and scalar screening masses are essentially 
degenerate.) The data show that the correlation length is dominantly driven by l/m,, much 
less by the gauge coupling /3. 
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FIG. 3: Correlation length (inverse pseudoscalar screening mass) vs inverse quark mass for Nj = 2 
sextet fermions on a 12^ x 6 volume. Bursts are (3 = 4.2; crosses, (3 = 4.4; squares, /? = 4.6, 
octagons, /? = 4.8, and diamonds, /? = 5.2. Lines connect points with the same /? values. 

B. Contrasting the deconfined phase for Nf = 2 fundamental fermions 



The reader might recall a fact about screening masses in the deconfined phase of funda- 
mental representation QCD: they scale with system size like 54 



H 



TT 



N.I 



(11) 



where is the lowest nonzero Matsubara frequency associated with antiperiodic boundary 
conditions in a lattice of temporal length Nt. The data in Fig. [3] are for Nf = 6, for which 
Eq. [TT] predicts mn = 27r/6 ~ 1 in the chiral limit. The limiting value of ttih is about half 
that value. In fact, Eq. [TT]does not reproduce any of the data. The situation at the one f3 
value (5.2) where I have many volumes, ranging from 12^ x 6 to 16^, is shown in Fig. HI The 
limiting correlation length does not seem to depend separately on A^^; rather, it appears to 
increase as the lattice volume increases. 



C. Scaling of correlation lengths? 

In the infinite volume limit we expect that the mass gap will disappear at vanishing 
fermion mass. The chiral susceptibility x = d'E{mq) / drriq is related to the volume integral 
of the isosinglet correlator C{x) = {'4)ip{x)il)ip{Qi)) . When the lightest mass in this channel 
vanishes, the susceptibility will diverge. As in QCD, this argues that the screening mass 
in this channel should vanish. Again, as in QCD, the calculation of isosinglet correlators 
is difficult due to the presence of disconnected diagrams, so as a first try I will mimic the 
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FIG. 4: Correlation length (inverse pseudoscalar screening mass) vs inverse quark mass for Nf = 2 
sextet fermions at /3 = 5.2. Diamonds, 12^ x 6 volume; octagons, 12^ x 8 volume, crosses, 12^ 
volume; squares, 16^ x 8 volume, and burst, 16^ volume. Lines connect points with the same P 
values. 

QCD analysis [s^ and look at the isotriplet scalar correlator. Just because the data is quite 
accurate, I will also look in the pseudoscalar channel. 

Fig. m shows us that data from different volumes seem to lie on a volume independent 
curve until the correlation length reaches a size comparable to the system size. Presumably 
the relation between rounding of the correlation length and volume is geometry dependent. 
So I will work with the two data sets I have, which have the same aspect ratio, 12^ x 6 and 
16^ X 8. I will try a fit of the form of Eq. [H varying the range and composition. 

A pair of representative fits are shown in Fig. O The left panel is the pseudoscalar and 
the right panel the scalar; I will define an exponent through l/rrij ~ (l/mg)^/^^ with j = Ti 
or oq. For the pictured fits the results are y.,, = 1.66(6), = 41.7/7 dof (pseudoscalar) 
and = 1.54(5), = 9.7/7 dof (scalar). It is relatively easy to find ranges of fits for the 
scalar mesons for which < 10 with l/yag drifting from .51(2) to .67(3) according to the 
fit range. The pseudoscalar data has smaller uncertainty and hence poorer confidence levels, 
X^/dof between 3-4 for 4-9 points, and 1/y ranging from 0.57(1) to 0.62(2). My data do not 
have the dynamic range to take my result too seriously. (If people had not spent the last 
30 years fitting ~ in QCD over the same range of masses, I would not have tried it.) 
Researchers with larger simulation volumes might like to attempt this, however. 
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FIG. 5: Representative fits to (a) pseudoscalar and (b) scalar correlation lengths. Data are dia- 
monds for 12^ X 6 volumes and squares for 16^ x 8 volumes. 

IV. DISCUSSION OF RESULTS FROM EIGENVALUES 
A. Volume scaling? 

Since I have argued that the gauge coupling appears to be irrelevant, results from any 
gauge coupling anywhere in the deconfined phase (in the basin of attraction of the FP?) 
should be equally meaningful. I chose a gauge coupling [3 = 5.2 and k = 0.1285. These are 
the second- smallest masses in Fig. HI At this gauge coupling simulations are easy to perform. 
At this hopping parameter the IR cutoff (lattice volume) dominates the effect of the quark 
mass (the AWI mass is about 0.042). I computed the values of the lowest 8 eigenvalues of 
the Hermitian squared Dirac operator D^D, which I then converted into eigenvalues of the 
overlap operator by the usual stereographic projection. All of the configurations collected 
at this parameter value have zero topological charge. Results were shown in Fig. [H 

Two pictures illustrate the quality of the data, combined into Fig [61 The left panel 
shows simulation time histories of the lowest four eigenvalues of the 16^ data set. Each 
measurement is separated by five HMC trajectories. The right panel shows the error on the 
average computed by blocking Nfy successive measurements together. The lowest eigenvalue 
is clearly the noisiest, but the autocorrelation time seems not to be too large: I bin two 
successive lattices {Nf, = 2 or At = 10 HMC units) together before averaging. 

The spectrum is nothing like the spectrum of eigenvalues from a theory with chiral sym- 
metry breaking, which consists of roughly equally spaced energy levels extending from the 
origin. The large spectral gap separating the lowest eigenvalue from the origin reflects the 
absence of chiral symmetry breaking in the weak coupling phase. 

I repeated the calculation of eigenvalues using fundamental-representation valence quarks 
from a fundamental-representation simulation in the deconfined phase. The gauge coupling 
is /3 = 5.5 and k = 0.126. This is a parameter set in which the screening masses have 
saturated at their Matsubara values. The picture is completely different - see Fig. [3 
Again, there is a spectral gap, but now it is dominantly affected by the short (antiperiodic) 
length of the lattice, Matsubara physics again. The eigenvalue distribution does not depend 
simply on L = V^^^. 
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FIG. 6: (a) Time history of overlap eigenvalues from the 16^ data set (in units of 5 HMC time 
steps), (b) Uncertainty on the average (Aj) as a function of bin size. Symbols are squares for i = 1, 
octagons for i = 2, diamonds for i = 3, and crosses for i = 4. 
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FIG. 7: Average value of ith eigenvalue of the fundamental-representation valence overlap operator, 
vs 1/L, where the lattice volume V is defined to be F = L^, from fundamental representation 
simulations. The actual volumes, moving from left to right across the graph, are 16^ x 6, 12^ x 8 
and 12^ x 6. The four lines are scaling curves {Xi)LP = constant, for p = 4, 3, 2, 1, from the top 
of the figure down. 

Now let us try to extract an exponent from Fig. [H This is not so easy: for theoretical 
input, all we have is Eq. [6l the finite-size scaling formula. One does not know a priori if 
it applies to all the eigenvalues or to only the lowest eigenvalues. One also does not know 
if there is some minimum volume for which it applies. (All of these questions have RMT 
answers for chirally-broken systems, but that is not what we are looking at here). So I will 
just proceed empirically: I will look at fits to individual eigenvalues, then groups of them. I 
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TABLE I: Exponent p from fits to individual eigenvalues, from the largest three volumes (16^, 
16^ X 8, 12"'), or the largest four (add 12^ x 8), or to all volumes (add 12^ x 6). 
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will fit all the data sets or drop smaller volumes and fit only the larger ones. 

I begin by fitting individual eigenvalues (lowest, second, and so on) to a power law, 
In(Aj) = Ai — pin L. I choose to fit to all five volumes, or four, or the largest three. The 
individual data points in each fit are uncorrelated, of course. Fits and chi-squareds are 
shown in Table [H Examples of fits are shown in Figs. [S] and M 

Clearly, I can also fit groups of eigenvalues. In that case, I fit In(Aj) = Ai — pin L for 
i = 1 . . . N eigenvalues. Now the data are correlated. I look at the quality of fits from 
uncorrelated fits, and then repeat by taking bootstrap averages of the data. The behavior 
is quite similar to the fits to individual eigenvalues. Some examples (with bootstrap errors) 
on the average, chi-squared from uncorrelated fits: 

• Fit the lowest 4 eigenvalues and biggest 3 volumes: p = 1.30(4), /dof = 25/(12 — 5) 

• Fit the lowest 4 eigenvalues biggest 4 volumes: p = 1.42(3), /dof = 67/ (16 — 5) 

• Fit all 8 eigenvalues and biggest 3 volumes: p = 1.46(3), x'^/dof = 331/ (24 — 9) 

• Fit all 8 eigenvalues and biggest 4 volumes: p = 1.43(2), x'^/dof = 372/(32 — 9) 

Examples of these fits are shown in Fig. [TUl 

Generally, the fits are poor but the trend of the four- volume fits is clear: p ~ 1.4. I 
cannot assign an error which is not a guess: there are too many unknowns. Are the volumes 
large enough? Are the data sets large enough? Nobody has tried such an analysis before. I 
think it is better to present my results and let others explore new directions. Nevertheless, 
if the reader desires a number with an error bar, p = 1.4(1) will capture the uncertainty 
over which fit to choose. 

B. Additional properties of eigenvalues 

Of course, there is more to the eigenvalue distribution than its average. In the absence of 
any theory, I will just compare cumulants, C{x) = n{x)/N where n{x) is the number of data 
points with a value smaller than x and N the total number of data points. I will rescale the 



13 



-In L 



-In L 




-2.3 ^ ' ... I .... I .... I ... . i 
-2.8 -2.7 -2.6 -2.5 -2.4 
-In L 




-2.8 -2.7 -2.6 -2.5 -2.4 
-In L 



FIG. 8: Three-volume fits to individual eigenvalues, (a) the lowest eigenvalue (b) the first excited 
state (c) the second excited state (d) the third excited state. The volumes are (from the left) 16^, 
16^ X 8 and 12'^. 

data at each volume by a factor [Li/LiY with p = 1.4 and Li = 16 in the plots. This is 
shown in FigfTTl 

Apart from the lowest eigenvalue the cumulants are very sharp, indicating a strongly 
localized eigenvalue spectrum. This was observed by the authors of Ref. [9|. The lowest 
eigenvalue has a broad distribution. This is not surprising; the spectral gap is large and the 
"pressure" on it from all the higher eigenvalues is asymmetric. Could the broad distribution 
be related to the relatively poor scaling fits for this mode? 

A direct calculation of p(A) involves histogramming the data. A possibly better technique 
is inspired by Ref. j56|, namely to integrate the density 

.A 

5(A) = / p{\)dX ~ A^+\ (12) 
Jo 

(basically the cumulant of all eigenvalues less than A). This was not a successful way to 
measure 7, since I did not have enough eigenmodes. Fig. [12] illustrates a typical data set, 
this time the 16^ x 8 one. A few minutes with a ruler will produce a power law in the vicinity 
of 3 for 1 < ^(A) < 8 (compared to 4/|/m = 2.9 for ym = 1.4) but this is clearly only of 
entertainment value: the curve flattens when the data sets are exhausted, and presumably 
at lower A values the higher (missing) modes would contribute as well. Ref. |56| analyzed 



14 




-S.8 -2.7 -2.6 -2.5 -2.4 -2.3 -2.8 -2.7 -2.6 -2.5 -2.4 -2.3 

-In L -In L 



FIG. 9: Four-volume fits to individual eigenvalues, (a) the lowest eigenvalue (b) the first excited 
state (c) the second excited state (d) the third excited state. The volumes are (from the left) 16^, 
16^ X 8, 12^, and 12^ x 8. 

far more eigenvalues. I can collect several tens of eigenvalues of the squared clover Dirac 
operator quite easily, but the analysis did not seem to be crisp due to the explicit chiral 
symmetry breaking in the action, shifting the zero, and I did not pursue this as far as I 
probably should have. 

At low A the plot shows a clear break. This is the part of the data which is completely 
dominated by the lowest eigenvalue, whose low-A edge is not given by a simple power law. 
There must be interesting physics here (as there is for eigenvalue distributions of chirally 
broken systems) but I do not know what it is. I also do not know if this is real physics or 
just statistics: if p(A) is small in some region of A, a big data set is needed to populate it. 

V. DISCUSSION 

Cautious readers might conclude that this paper is content-free: a curious Fig. [T] which 
appears only at the end of a long set of improbable constructions. In fact, I present the 
paper because Fig. [T]is so striking. Most of my analysis assumed that the underlying theory 
was conformal in the massless, infinite volume limit. I believe this assumption but it needs 
verification. Whether it is or not, one fact seems true: The deconfined phase of SU{3) gauge 
theory with Nf = 2 flavors of sextet fermions seems to have quite different properties from 
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FIG. 10: Combined fits to several eigenvalues, (a) Three-volume fits to the 8 lowest eigenvalues, 
(b) four volume fits to the lowest 8 eigenvalues, (c) Three-volume fits to the lowest 4 eigenvalues, 
(d) four volume fits to the lowest 4 eigenvalues. The volumes are (from the left) 16^, 16^ x 8, 12^ 
and in (b) and (d) 12^ x 8. 



the deconfined phase of SU{3) gauge theory with Nf = 2 flavors of fundamental fermions. 

Let us assume that the exponent I was measuring is truly what I claim it is: the leading 
exponent of a theory with one relevant coupling (the mass) which is critical at zero mass. As 
far as a lattice simulator can tell from reading the continuum literature, the answer i/^ ~ 1-4 
is very interesting: scaling exponents far from their engineering dimensions seem to be much 
sought after. (A partial subset of the literature is Refs. 57|, [SSl, l59|, [GO, l6l|, l62(].) Recalling 
that the exponent of the condensate is 7 = 4/?/m — 1, this is an exponent of about 7 = 1.9(2). 

A common prediction in the literature is that a lower bound on 7 is 7 = 2 or j/^ = 4/3. 



The earliest appearance of this result I can find is by Cohen and Georgi [57[, who obtain it 



by solving a gap equation. In a recent paper Kaplan, Lee and Son argue on different 
grounds that 7 should always exceed 2. My result is barely consistent with this prediction, 
with an uncertainty which in the end seems a bit inadequate. 

Ryttov and Sannino ^] have a supersymmetric QCD - inspired beta function for gauge 
theories with higher dimensional representations of fermions. The anomalous dimension is 
predicted to be 

11C2(G) -AT(R)Nf 13 , , 

7 = 3 ' ^ = 3 - — = 1.7. (13) 
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FIG. 11: Cumulants, where the eigenvalue is rescaled by {Li/LiY with p = 1.4 for the four largest 
volumes: (a) 16^ (no rescaling) (b) 16^ x 8 (c) 12^ (d) 12^ x 8. 
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FIG. 12: ^(A) vs A (defined in Eq. ^ from the 16^ x 8 data set. 
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There are some models of technicolor where the new physics sector is conformal or nearly 
so (for example [g^, Is^); large anomalous dimensions (7 near 1) are desirable to achieve 
their phenomenological goals. This does not seem to be the case for sextet QCD. 

The theoretical literature almost always addresses many models simultaneously, and val- 
idating (or otherwise) their assumptions requires data from many systems. ?/m's for other 
lattice IRFP theories are needed. 

From a simulation point of view, the present study could be improved by eliminating 
partial quenching. Perhaps the best candidate to do this is SU (2) color with adjoint fermions. 
Codes may be fast enough to allow for simulations with dynamical overlap fermions on 
reasonable volumes. 

A limitation of the use of eigenvalues is that at present, the theory is just finite size scaling. 
RMT predictions were rich because the number of theoretically motivated observables was 
large. Analysis like what I am doing needs additional theoretical input. 

The published simulations of Nf = 12 fundamental fermions use staggered fermions. 
The technology for measuring susceptibilities for the condensate in the context of QCD 
thermodynamics is quite well developed. It would be interesting to apply it to present 
needs. Staggered fermion eigenvalues are also cheap to compute but explicit flavor breaking 
effects could hamper their interpretation. 

Ultimately, the only universal quantities associated with IRFP theories are their critical 
exponents. Can lattice simulations measure them? 
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